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The quasiparticle model is extended to investigate the properties of strange quark matter in a strong magnetic 
field at finite densities. For the density-dependent quark mass, self-consistent thermodynamic treatment is 
obtained with an additional effective bag parameter, which depends not only on the density but also on the 
magnetic field strength. The magnetic field makes strange quark matter more stable energetically when the 
magnetic field strength is less than a critical value of the order 10^ Gauss depending on the QCD scale A. 
Instead of being a monotonic function of the density for the QCD scale parameter A > 126 MeV, the effective 
bag function has a maximum near 0.3 ^ 0.4 fm^-'. The influence of the magnetic field and the QCD scale 
parameter on the stiffness of the equation of state of the magnetized strange quark matter and the possible 
maximum mass of strange stars are discussed. 

PACS numbers: 24.85.+p, 12.38.Mh, 21.65.Qr, 25.75.-q 



I. INTRODUCTION 



Since strange quark matter (SQM) was speculated by Witten as the possible true ground state of strong interaction matter 
lit], the properties of SQM in bulk, as well as in finite size, the so called strangelets, have been extensively studied in thepast 
decades JJ-qII. The new form of matter is possibly produced by terrestrial relativistic heavy-ion collision experiments fl] or 
exists in the interior of compact stars It was found that the stability of SQM is strongly affected in a strong magnetic field 
Ht]. The large magnetic fields in nature are normally associated with astrophysical objects, where the density is much higher 
than the nuclear saturation. The typical strength could be of the order ^ lO'^ G on the surface of pulsars i^. Some magnetars 
can have even larger magnetic fields, reaching the surface value as large as lO''* ~ lO'^ G ifioll . In the interior of compact stars, 
the maximum possible magnetic field strength is estimated as high as ~ lO'^ G. The origin of the strong magnetic fields can be 
understood in two ways. One is the amplification of the relatively small magnetic field during the star's collapse with magnetic 
flux conservation ilUl . The other is the magnetohydrodynamic dynamo mechanism with large magnetic fields generated by 
rotating plasma of a protoneutron star llT2ll . 

Because a strong magnetic field influences the single particle spectrum while all quarks are charged, SQM in the inner part of 
a compact star may show specific properties. Specially, for example, the strong magnetic field leads to a more stable polarized 
strange quark star (SQS)|l3l]. In heavy-ion collisions experiments, the magnitude of a magnetic field plays an important role in 
studying the deconfinement and chiral phase transitions. In the LHC/CERN energy, it is possible to produce a field as large as 
5 X 10''* G H. 

With various phenomenological confinement models, many works on the properties of magnetized SQM(MSQM) have been 
done by a lot of researchers. Based on the conventional MIT bag model, quark matter in a strong magnetic field was studied 
by Chakrabarty [j^, and significant effect on the equation of state had been found. Furthermore, the magnetized strangelets at 
finite temperature was investigated by Felipe et. al. in their recent work lfl5l[l6ll . In Ref. iH, the effect of an external magnetic 
field on the chiral dynamics and confining properties of SQM were discussed in the linear sigma model coupled to the Polyakov 
loops. The special properties of MSQM were also investigated with the Nambu-Jona-Lasinio (NJL) model lfl8l- l2l]| . The MIT 
bag model, the two-flavor NJL model, and the chiral sigma model had also been compared in studying the MSQM ll22ll . 

In literature, the quasiparticle model, where the effective quark mass varies with environment, was also successfully employed 
by many authors to study the dense strange quark matter in the absence of an external magnetic field ll23l - l25ll . The main advantage 
of the quasiparticle model is that it can explicitly describe quark confinement and vacuum energy density for bulk matter ll24ll and 
strangelets 126(1 . The aim of this article is to extend the quark quasiparticle model to studying the magnetized quark matter. We 
find a density- and magnetic-field-dependent bag function. Accordingly, a self-consistent thermodynamic treatment is obtained 
with the new version of the bag function. The effect of a magnetic field on the bag function and the stability of MSQM will be 
discussed. It is found that the magnetic field makes SQM more stable when the magnetic field strength is less than a critical 
value of the order 10^ G depending on the QCD scale A. 
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This paper is organized as follows. In Sec. 2, we derive the thermodynamic formulas in the quasiparticle model when the 
magnetic field becomes rather important, and then demonstrate the effective bag function for the case of both constant and 
running coupling, respectively. In Sec. 3, the stability properties of MSQM, the effective bag function, and the mass-radius 
relation of magnetized quark stars are investigated, and discussions are shown about the effect of the magnetic field and QCD 
scale parameter. The last section is a short summary. 



II. THERMODYNAMIC TREATMENT IN A STRONG MAGNETIC FIELD 




The important feature of the quasiparticle model is the medium dependence of quark masses in describing QCD nonperturba- 
tive properties. The quasiparticle quark mass is derived at the zero-momentum limit of the dispersion relations from an effective 
quark propagator by resuming one-loop self-energy diagrams in the hard dense loop (HDL) approximation. In this paper, the 
effective quark mass is adopted as ll24i uTI l28ll 

'«'(mO-^ + \/^ + ^, (1) 

where m,o and /i, are, respectively, the quark current mass and chemical potential of the quark flavor /. The constant g is the 
strong interaction coupling. One can also use a running coupling constant g{Q/A) in the equations of state of strange matter 
instead of a constant g 12911 . In our recent work by using phenomenological running coupling [|2SI, the quark masses were 
demonstrated to decrease with increasing densities at a proper region. 

Here, we assume the g value is in the range of (0,0.5), as done in the previous work i24ll . The current mass can be neglected 
for up and down quarks, while the strange quark current mass is taken to be 120 MeV in the present calculations. Because the 
vanishing current mass is assumed for up and down quarks, Eq. ([T]) is reduced to the simple form 

Instead of inserting the effective mass m, directly into the Fermi gas expression, we will derive the expressions from the self- 
consistency requirement of thermodynamics. The quasiparticle contribution of the flavor / to the total thermodynamic potential 
density can be written as 



diT 



jJ^°°{ln[l+e-(^''.-M,)/r 

l+e-fep+M,)/r]|d3^;j, (3) 



where T is the system temperature and c/, is the degeneracy factor [di — 3 (color) for quarks and di — 1 for electrons]. All the 
thermodynamic quantities can be derived from the characteristic function by obeying the self-consistent relation ll30ll . 

To definitely describe the magnetic field of a compact star, we assume a constant magnetic field (B,„,^ = B,„) along the z 
axis. Due to the quantization of orbital motion of charged particles in the presence of a strong magnetic field, known as Landau 
diamagnetism, the single particle energy spectrum is IBlll 



Ei ^\IpI + mj + eiB,„ {In +.?+!), (4) 

where is the component of particle momentum along the direction of the magnetic field e,- is the absolute value of the 
electronic charge (e.g., e,- = 2/3 for the up quark and 1/3 for the down and strange quarks), n ~ 0, 1,2,..., are the principal 
quantum numbers for the allowed Landau levels, and s = ± 1 refers to quark spin-up and -down state, respectively. For the sake 
of convenience, we set 2v = 2« + s + 1, where v = 0, 1,2, .... The single particle energy then becomes [Ql 



Ei = \Jpl + m} + 2veiB,n. (5) 
On application of the quantized energy levels, the integration o\ei dpxdpy in Eq. (O is replaced by the rule 

dpxdpy^2neiB,„ ^ (6) 

s=±l n 
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Because there is the single degenerate state for v = and the double degenerate state for v 7^ 0, we assign the spin degeneracy 
factor (2 — Syo) to the index v Landau level. The thermodynamic potential density of Eq.© in the presence of a strong field can 
thus be written as 



n,(r,m„MO - -7-^ 1(2- 5vo) {ln[l +exp(^)] +ln[l +exp(-^)]}rfp. 

At zero temperature, Eq. Q is simplified to give 

deB r\/Jif-MW 



(7) 



(8) 



where My = y 'wf + 2 Ve,B„, is the quark effective mass in the presence of a magnetic field. In the case of zero temperature, 
the upper limit Vmax of the summation index v can be understood from the positive value requirement on the logarithm and 
square-root function in Eq. dU. So we have 

V<V,„„x = mt[ ^' ' ], (9) 

where "int" means the number before the decimal point. 

Accordingly, the pressure P, the energy density E, and the free energy density F for SQM at zero temperature read 13 



P = -Q.f-B*, (10) 
E = F ^Q.f + Y,n,ni + B*. (11) 

i 

Here Q.f = is the free quasiparticle contribution with the summation index going over all flavors considered. The notation 
B* denotes the effective bag function and it can be divided into two parts: jU,-dependent part and the definite integral constant 
part, i.e., B* = +Bo (' = u, d, and s) where Z?o is similar to the conventional bag constant and Bi{pii) is the chemical 

potential dependent function to be determined. 

The derivative of the thermodynamic potential density H, with respect to the quark effective mass m, has an analytical expres- 
sion, i.e.. 

The quark particle number density of the component / is given as 

dieiB, 



^(2-5vo)VA''-< (13) 



'^'^ v=0 

In the literature, there are three methods to construct a consistent set of thermodynamical functions with the effective quark 
masses. One is applied in the quark mass density-dependent model in Refs. lf33l[34ll . where all thermodynamic quantities are 
derived by direct explicit function and implicit function dependent relations. The second is the treatment in the NJL model, 
where the dynamical quark masses are solutions of the gap equation coupling the quark condensates ll2ll[35ll . The energy and 
pressure functions are modified accordingly. The third method is to get a self-consistent thermodynamical treatment with an 
effective bag constant to describe the residual interaction ||36ll. The effective bag constant acts as a part of a modified pressure 
function. Here, we employ the third method. The following requirement is introduced and applied as in Refs. ll24il37ll . 



dP 
dnii 



0. (14) 



From a physical viewpoint, the constraint can make the formula of particle number function consistent with standard statistical 
mechanics. From Eqs. ( fTOl i and (fTTT i. it can be understood that the effective bag constant leads an additional term in the 
modification in the energy and pressure functions. 
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Considering Eq.(fl4]i. we have the vacuum energy density Bi{fXi) through the following differential equation 

dBi(^d^^_5% 
d/Z; dm, dnii 

If we assume the vanishing current quark mass, one can integrate Eq. ( fTSl l under the condition = 0) = and have 



(15) 



dnii 



27r2 



v=0 



'■)2 



(0 



-)d^i, 



where the lower limit of the integration over jXi is different from that in Ref. Il24tl . Its critical value /if should satisfy 

jUf2-m^-2ve,B„, >0 



(16) 



(17) 



To reflect the asymptotic freedom of QCD, the calculation must be changed by including the running coupling constant. The 
approximate expression for the running quantity g{n) reads Issll . 



^'(r = o,M) 



48 TT^ 

~29' 



ln( 



0.8/12 



-1 



(18) 



where A is the QCD scale parameter, the only free parameter in the theory determined by experiments. The magnitude of A 
controls the rate at which QCD coupling constant runs as a function of exchanged momentum (see Ref. 1221). After applying 
the running coupling constant ( fTSl l, the effective bag function in Eq. ( fT6b is changed into 



B^{}^i. 



do. 



C'W; 



dm, 
d/ii 



deB f^' 



v=0 



(0 



dju. 



(19) 



where the lower limit of the integration /if satisfies Bi{jj.^ = 0). Differently from the constant coupling case, the critical value 
/if can be obtained by inserting the running coupling constant in Eq. ( fTSl l into the condition ( fTTI ). The value of /if depends not 
only on the chemical potential of quarks but also on the Landau energy level. 



III. PROPERTIES OF MAGNETIZED STRANGE QUARK MATTER 

In this section, the properties of MSQM are studied with the new version of the quasiparticle model in the presence of a strong 
magnetic field. We will investigate the properties with a density- and magnetic-field-dependent bag function. Then, we discuss 
the effect of the QCD scale parameter and the strong magnetic field on the effective bag function and strange quark stars. 



A. The stability property of bulk magnetized SQM 

As is usually done, the SQM is treated as a mixture of u-, d-, s- quarks and electrons with neutrinos entering and leaving the 
system freely. To obtain the equations of state (EoS) of magnetized SQM, a set of equilibrium conditions-the weak equilibrium, 
baryon number conservation, and electric charge neutrality-should be considered by the following relations [HI l39l - l4Tll : 

/i„ + /ig = /i,/ = /lj, (20) 
«M + nrf + «.s = 3nB, (21) 

:^nu - -n,i - -n, - = 0. (22) 

Equation ( |20b is the chemical equilibrium condition maintained by the weak-interaction processes such ass + u—>-u + d 
and s ^ u + e + Ve etc., Eq. i2l[ is from the definition of the baryon number density «b, and Eq. ( |22] | is the charge neutrality 
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condition. For a given baryon number density «b, we can obtain the four chemical potentials /i„, ji^, ji^, and /Ze by solving the 
four equations in Eq. (|20li-(l22li. Other thermodynamic quantities, such as the energy density and pressure, can then be calculated 
from the formulae derived in the previous section. A little difference is that the Maxwell contribution has been included in our 
numerical calculations, i.e., the quasiparticle contribution Q.f is replaced by ll42l - l44ll 

n = % + (23) 

where the second term is the pure Maxwell contribution of the magnetic field itself. 

In Fig.[Tl the energy per baryon of MSQM is shown as functions of the density for several g values. For comparison purposes, 
we have also plotted the previous results in Ref. ll26ll by setting B,„ = 0. The solid curves are for MSQM, while the dotted ones 
are for the corresponding nonmagnetized SQM. The two groups of curves have apparently similar density behavior. Obviously, 
however, the MSQM has lower energies than the nonmagnetized SQM. To show the effect of different coupling constants, we 
adopt three values of g. 
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FIG. 1: The energy per baryon versus the density at fixed coupling constant g = 2,3,4 for magnetic field strength B,„ = lO'^G. Compared 
with the nonmagnetized strange quark matter (the dotted curves with B,„ = 0), the magnetized case has a lower energy per baryon. 

In the quasiparticle model, the parameter g stands for the coupling strength and it is related to the strong interaction coupling 
constant aj hy g = y/Ana's. Therefore, the g value has a large effect on the stability of SQM [|45|| . To satisfy the requirement of 
QCD asymptotic freedom, the running property of the coupling parametrization should be considered. In Fig. |2] we show the 
running coupling constant as functions of the baryon number density ns- The three lines are obtained with different values of 
A. It is very obvious from Fig. |2]that the running coupling g is a decreasing function of the density. With a bigger A value, the 
coupling g is also bigger at any fixed density. 

In Fig. [3] we show the same quantities as in Fig, [T] with the running coupling constant, respectively for the two values of the 
different magnetic field 10^^ G (dashed lines) and lO'^G (solid lines). It is clearly seen that the energy per baryon increases with 
increasing the QCD scale parameter A, i.e. SQM has a lower energy per baryon with smaller A value at a fixed strong magnetic 
field. This effect of the QCD scale parameter is consistent with the constant coupling case in Fig. [T] because larger A means 
bigger coupling as indicated by Eq. ( fTSl l. 

An obvious observation from Fig.[3]is that there is a minimum energy per baryon for each pair of the parameters A and B„,. In 
Fig.m therefore, we show how the minimum energy of MSQM varies with the magnetic field strength. The QCD scale parameter 
is taken to be 180 MeV (the upper dashed curve) and 120 MeV (the lower solid curve) respectively. It is found on each curve 
that there is another minimum value corresponding to a critical magnetic field strength For the values of A = 120 MeV and 
180 MeV, the corresponding BJ„ equals 2.15 x lO'^ G and 2.34 x lO'^ G respectively. When the magnetic field strength is less 
than fiJJ,, the minimum energy per baryon decreases with increasing the strength of the magnetic field. When the magnetic field 
strength exceeds or, equivalently, when the magnetic energy scale approaches the QCD scale, i.e., y/eB,„ ~ 76.9 MeV, the 
field energy itself will have a considerable contribution to the energy of SQM and hence the energy per baryon increases with 
the magnetic field strength. In Fig. [3] the magnetic field strength is taken to be the corresponding critical value. 

Because we study magnetized strange quark matter in the "unpolarized" approximation, it is appropriate to estimate the 
maximum magnetic field strength when such an approximation can be reliable. To this end, in principle, we can investigate the 
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FIG. 2: The running coupling constant g versus tlie baryon number density at different A values with the magnetic field B,„ = lO'^G. The 
upper lines correspond to larger values of A. 
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FIG. 3: The energy per baryon £/;7b of stable MSQM versus the number density at different A values with the corresponding critical magnetic 
field strength BJJ,. 



polarized quarks with spin up (+) and down (-) by introducing the polarization parameter ^; as ifisi l4lll 



"h - (24) 



where n.^' and n[ ' denote the number density of spin-up and -down /-type quarks. For the sake of simplicity, we assume a 
common polarization rate ^ for u-, d-, and i-quarks, i.e., = = = ^- In SecUll the summation for fixed spin s = +1 or 
s = — 1 should go over the principal quantum numbers n instead of v. The degeneracy factor (2 — 5vo) in Eqs.d?), (O, (fT2] i and 
(T3[ should be deleted because the spin degeneracy disappears for polarized particles. The polarization parameter < ^ < 1 
will decrease with increasing the number density. Assuming a larger value of the polarization £, = 0.6, the energy is enlarged by 
4.5%. In fact, even for v ery larger magnetic field B,„ = 5 x 10'**G , the parameter ^ remains in the range (0.01 ~ 0.02) when the 
density ng > 0.2 fm llSll . We do the numerical calculation and find that the free energy per baryon will be enlarged by 0.8% 
at ^ = 0.1. So the effect of the unpolarized approximation on the discussion of the stability of SQM is very small, especially 
when the magnetic strength is less than lO^^G which is an estimated maximum possible strength of the interior magnetic field. 
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FIG. 4: The energy of magnetized strange quark matter varies with the magnitude of the strong magnetic field for the fixed QCD scale 
parameter A = 180 MeV (the upper dashed curve) and A = 120 MeV (the lower solid curve) respectively. With decreasing the magnetic field 
strength, the energy per baryon approaches gradually to the value without a magnetic field indicated by a horizontal dashed-dotted line. 



B. The effective bag function for magnetized SQM 

The effective bag function B* is generally used to represent the vacuum energy density for dense QCD matter ll46ll . Com- 
paring it with the standard statistical mechanics, one can recover the thermodynamics consistency of system density- and/or 
temperature-dependent Hamiltonian with the extra term B* . The meaning of B* plays an important role in studying properties of 
quark matter. The interpretation of B* was first given by Gorenstein and Yang in Ref . IjtIi . In quasiparticle model, because the 
dispersion relation is density- and/or temperature-dependent, B* is regarded as the system energy in the absence of quasiparticle 
excitations, which cannot be discarded from the energy spectrum 14711 . In this sense, B* acts as the bag energy or bag pressure 
through the application in bag-like model. One can interpret the confinement mechanism considering B* as the difference of 
perturbative vacuum and physical vacuum. 

In addition to the constant value Bq of the bag model, the expression of B* has been developed in several different forms. Li, 
Bhalerao, and Bhaduri obtained the temperature-dependent bag constant in the QCD sum-rule method ll48ll . Song obtained a 
/i- and T-dependent bag constant by incorporating one-loop correction in imaginary time formulation of finite temperature field 
theory IH, 

^* ^) = ^0 - tl^^' + + • (25) 

In the work of Burgio ifsoll . the Gaussian parametrization of density dependence of B* is employed as, 

=Boo+(Bo-Bo„)exp(-7(«B/«of ), (26) 

where the parameters Boa, y, and «o are given in Ref. ifsoll . The effective bag constants in these previous works are all monoton- 
ically decreasing functions of the density and temperature ifsill . In our present work, the effective bag function B* is associated 
with a magnetic field and consequently has a different density behavior. We thus plot the effective bag function B* versus the 
baryon number density with different A values in Fig.|5] The dashed lines are for the magnetic field strength Z?,„ = lO'^G, while 
the solid lines are for a higher magnetic strength B^ = lO'^G. The open circles indicate nonmagnetized SQM. The numerical 
results show an important property that the effective bag function B* remains decreasing monotonously with increasing densities 
for smaller A = 120 MeV. But for larger value A = 180 or 200 MeV, the bag function B* has a maximum value at about 2^3 
times the nuclear saturation density 0.16 fm^^. Generally, when the QCD scale parameter is bigger than the critical value 126 
MeV, the effective bag function is not a monotonic function and reaches a maximum value BJ^^x the density range 0.3 ^ 0.4 
fm-3. 

Since the QCD scale parameter A plays a great role on the effective bag function B*, we plot the bag function B* of stable 
SQM, i.e., P = 0, versus A on the left axis in Fig,|6l If one requires that the bag function B* should be a nonmonotonic decreasing 
function of the density, the A value should be bigger than the critical value 126 MeV. The corresponding baryon number density 
iiB marked by a dashed line on the right axis is also plotted. The bag function B* and the baryon number density ns all increase 
with the QCD parameter A. 



8 




FIG. 5: The effective bag function B* for MSQM versus tlie baryon number density at different A values. B„, = lO'^G and B„, = lO'^G are 
marlced by dashed lines and solid lines, respectively. Only for larger A, the B* has a maximum value SJi^x indicated by asterisks at 2 ~ 3 times 
nuclear saturation density. 
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FIG. 6: A dependence of the bag function B* (solid line) and baryon number density (dashed line) corresponding to the zero pressure. The A 
should be larger than the critical value 126 MeV to produce a nonmonotonic behavior of B*. 



C. Mass-radius relation of magnetized strange quark stars 



Strange quark stars, a family of compact stars consisting completely of deconfined «, d, s quarks, have attracted a lot of 
researchers. The gravitational mass (M) and radius (R) of compact stars are of special interest in astrophysics. The strange 
quark stars were studied by many authors as self-bound stars different from neutron stars. It is pointed out that the possible 
configuration of compact stars, such as the strange hadrons, hyperonic matter, and quark matter core, can soften the equations 
of state of neutron stars lf52l - l54ll . In this section, we calculate the mass-radius relation of magnetized SQS together with the 
effective quark mass scale. Using the EoS of MSQM in the proceeding sections, we can obtain M and R by numerically solving 
ToIman-Oppenheimer-Volkoff(TOV) equations when fixing a central pressure P^. Varying continuously the central pressure, we 
can obtain a mass-radius relation M(R) in Fig. [T] The stable branches of the curves must satisfy the condition dM/dPc > 0. In 
this way, we can find the maximum mass along the same curve, which is denoted by full dots in Fig.|7] Other solutions, on the 
left side of the maximum mass, are unstable and collapsible. 

It is seen from Fig. |2]that the maximum mass is bigger with a smaller A value and an extremely large magnetic field. However, 
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it is still not as big as the recently observed maximum mass of PSR J1614-2230 llssll. This may mean that a simple ordinary 
phase cannot explain the large mass. Some new phases, e.g., the superconductivity phase in dense matter 1156145811 . should be 
further studied in the future. 
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FIG. 7: The mass-radius relation of SQS at different A values with different magnetic fields B„, — lO'^G (dashed lines) and B,„ = lO'^'O (solid 
lines). The maximum masses on all curves are marked by full dots. 



IV. SUMMARY 

We have extended the quark quasiparticle model to study the properties of strange quark matter in a strong magnetic field at 
finite density. The self-consistent thermodynamic treatment is obtained through an additional bag function. The bag function 
depends not only on the quark chemical potentials but also on the magnetic field strength B„,. By comparison with the nonmag- 
netized quark matter, we find that the magnetic field can enhance the stability of SQM when the magnetic field strength is lower 
than a critical value of the order lO'^ G. But when the magnitude of the magnetic field is larger than the critical value the 
magnetic energy will have a considerable contribution to the energy of SQM. So the energy per baryon of MSQM increases with 
increasing the field strength. Because the quark masses depend on the corresponding chemical potential, an additional effective 
bag function, which depends not only on the chemical potentials but also on the magnetic field strength, appears in both the 
energy density and pressure. The effective bag function has a maximum at about 2^3 times the saturation density when the 
QCD scale parameter is larger than 126 MeV. Although an unpolarized approximation is assumed, we find the energy per baryon 
would increase by 0.8% for the usual polarization parameter when riB > 0.2 fm^^. 

On application of the new equation of state of the magnetized strange quark matter in ordinary phase to calculate the mass- 
radius relation of a quark star, it is found that the maximum mass does not explain the the newly observed maximum mass of 
about two times the solar mass. This means that other phases, e.g. superconductivity and/or mixed phases, might be necessary 
to explain the new astronomic observations, and further studies are needed. 
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